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ABSTRACT 

The luminous infrared galaxy Arp 299 (IC 694+NGC 3690) is studied using optical integral field 
spectroscopy obtained with the INTEGRAL system, together with archival Hubble Space Telescope 
WFPC2 and NICMOS images. The stellar and ionized gas morphology shows A-dependent variations 
due to the combined effects of the dust internal extinction, and the nature and spatial distribution 
of the different ionizing sources. The two-dimensional ionization maps have revealed an off-nuclear 
conical structure of about 4 kpc in length, characterized by high excitation conditions and a radial 
gradient in the gas electron density. The apex of this structure coincides with Bl region of NGC 3690 
which, in turn, presents Seyfert-like ionization, high extinction and a high velocity dispersion. These 
results strongly support the hypothesis that BI is the true nucleus of NGC 3690, where an AGN 
is located. In the circumnuclear regions Hii-like ionization dominates, while LINER-like ionization 
is found elsewhere. The Ha emitting sources with ages from 3.3 to 7.2x10^ years, have masses 
of between 6 and 680x10^ M© and contribute (extinction corrected) about 45% to the bolometric 
luminosity. The ionized (Ha) and neutral (NaD) gas velocity fields show similar structure on scales of 
several hundred to about 1 kpc, indicating that these gas components are kinematically coupled. The 
kinematic structure is complex and on scales of about 0.2 kpc does not appear to be dominated by the 
presence of ordered, rotational motions. The large velocity dispersion measured in NGC 3690 indicates 
that this galaxy is the most massive of the system. The low velocity amplitude and dispersion of the 
interface suggest that the ionized gas is slowly rotating or in a close to quiescent phase. 
Subject headings: galaxies: individual (Arp 299, IC 694, NGC 3690) — galaxies: interaction — 

galaxies: starburst — galaxies: nuclei — galaxies: stellar content — galaxies: 

kinematics and dynamics — infrared: galaxies 



1. INTRODUCTION 

Arp 299 is a nearby luminous infrared galaxy 
(d-43 Mpc for Ho=70 km s-^Mpc"!, 0.7, 
flA— 0.3) with an infrared (8-1000 fim) luminos- 
ity of L/fl=5.7x lO^^Lo^. It is an interacting sys- 
tem (Fig. ^ composed of two individual galaxies 
(IC 694-KNGC 3690"^) in an early dynamical stage. Hi- 
bbard & Yun (1999) have discovered an Hi tidal tail 
of 180 kpc in length, with no evidence for tidal dwarf 
galaxies. Powerful starburst regions with star formation 

^ Based on observations with the William Herschel Telescope 
operated on the island of La Palma by the ING in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de As- 
trofisica de Canarias. Based also on observations with the NASA- 
ESA Hubble Space Telescope, obtained at the Space Telescope and 
Science Institute, which is operated by the Association of Universi- 
ties for Research in Astronomy, Inc. under NASA contract number 
NAS5-26555. 

^ Departamento de Astroffsica Molecular e Infrarroja, Instituto 
de Estructura de la Materia, CSIC, E-28006 Madrid, Spain. 

^ Space Telescope Science Institute, 3700 San Martin Drive, Bal- 
timore, M.D. 21218, USA. 

On leave from the Instituto de Astrofisica de Canarias (lAC) 
and from the Consejo Superior de Investigaciones Cienti'ficas 
(CSIC), Spain. 

Instituto de Astroffsica de Canarias (lAC), 38205 La Laguna, 
Canary Islands, Spain. 

^ IRAS fluxes obtained from Moshir et al. (1993). L//j derived 
from the equation detailed on Sanders & Mirabel (1996). 

Through this paper, we use the notation introduced by Gehrz 
et al (1983) for the two components of Arp 299: the nucleus of 
IC 694 (eastern component) is region A, while the main sources in 
NGC 3690 (western component) are Bl, B2, C(C1-|-C2-|-C3) and 
C (see also Alonso-Herrero et al 2000). 



rates of up to about 100 Mq yr"^ have been identified 
(Alonso-Herrero et al. 2000, hereafter AAHOO). More- 
over, the system presents several regions where supernova 
remnants have been detected at radio wavelengths (Neff 
et al. 2004). More recently, Mattila et al. (2005) have 
discovered a new supernova using near-IR images, rein- 
forcing the prediction that this galaxy has a high rate of 
supernova explosions (Mannucci et al. 2003). 

Optical spectroscopic studies have classified IC 694 as 
starburst, and NGC 3690 as starburst/LINER (e.g., Co- 
ziol et al. 1998), and mid-IR observations classified the 
system as starburst (e.g., Laurent et al. 2000). Using 
Hubble Space Telescope (HST) imaging and near-infrared 
spectroscopy, AAHOO modeled the star formation prop- 
erties of the system without the need for an AGN. The 
first hint for the presence of an AGN in Arp 299 came 
from X-ray observations with the BeppoSAX satellite by 
Delia Ceca et al. (2002). These authors concluded that 
a highly absorbed (N^ ~ 2.5x10^** cm"^) AGN, with an 
intrinsic luminosity of Lo.5_ioofeey — 1.9x10^"^ erg s"^, 
is present in this galaxy, although due to the low spatial 
resolution its exact location could not be determined. 
Using Chandra and XMM-Newton observations, recent 
studies (Ballo et al. 2004) have inferred the presence of 
two AGNs, one located in NGC 3690 (source Bl; 6.4 keV 
line), and another in IC 694 (source A; 6.7 keV Fe-Ha 
line). In addition, Gallais et al. (2004) suggest the pres- 
ence of an AGN in NGC 3690 based on the shape of the 
mid-IR continuum, while in IC 694 they only detect spec- 
tral features characteristic of massive, dust enshrouded 
star formation. 
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This paper presents a detailed study of the two- 
dimensional star formation, kinematics, and ionized gas 
properties of Arp 299, covering the two individual galax- 
ies of the system and the interface area. At a distance of 
about 43 Mpc, Arp 299 is one of the nearest interacting, 
luminous infrared galaxies, and therefore these observa- 
tions provide high spatial resolution (about 0.205 kpc 
arcsec"^) and allow us to investigate in detail the struc- 
ture and physical properties of this galaxy. The paper is 
organized as follows: §2 describes the observations and 
data reduction, while the selection of regions under study 
is described in §3. The next sections present a discussion 
of the results, with special emphasis on the stellar and 
ionized gas morphology O, the internal dust distribu- 
tion and extinction effects (311), the extinction estimates 
using optical and near-IR emission lines f ii5.4|) . the two 
dimensional ionization structure (Ijnjl, the electron den- 
sity (33), the origin of the LINER-hke ionization (iQ, 
the star- forming properties O, and the kinematics and 
dynamical mass derivations f tllO() . Finally ill II presents 
a short summary of the main new results. 

2. OBSERVATIONS AND DATA REDUCTION 

Integral Field Spectroscopy (IFS) of Arp 299 was ob- 
tained on 2004 January 15 and 17 using INTEGRAL, 
a fiber-based integral field system (Arribas et al. 1998) 
connected to the Wide-Field Fiber Optic Spectrograph 
(WYFFOS) (Bingham et al. 1994) and mounted on the 
4.2 m William Herschel Telescope. Two different INTE- 
GRAL configurations were used for observing the sys- 
tem (see pointings in Fig. The fiber bundle with 
the largest field of view (FoV), SB3, was centered on 
a position equidistant from the centers of the galaxies 
and therefore used to map the diffuse, interface region 
between the two galaxies (IC 694-HNGC 3690) that form 
the Arp 299 system, as well as the galaxy external regions 
at distances of about 4 kpc. The SB3 bundle consists of 
115 fibers each 2'.'7 in diameter, covering a rectangular 
area of 33'.'6x29'.'4 on the sky. In addition 20 fibers, 
forming a circle of 45'.'0 in radius, were used to obtain 
a spectrum of the sky simultaneously. The SB2 bundle 
was used to map each of the two galaxies (IC 694 and 
NGC 3690, eastern and western galaxies, respectively on 
Fig. ^ individually. It consists of 189 fibers each 0'.'9 in 
diameter, covering a rectangular area of 16'.'0x 12'.'3 on 
the sky. In addition 30 fibers, forming a circle of 45'.'0 in 
radius, were used to measure the sky simultaneously. 

The spectra were taken using a 600 lines mm^^ grating, 
with an effective spectral resolution (FWHM) of approx- 
imately 6.0 A and 9.8 A for the SB2 and SB3 bundles re- 
spectively, and covering the spectral range AA4400— 7000 
A. The total integration time was 3600 s for each point- 
ing with the SB2 bundle, split into three individual expo- 
sures of 1200 s each. The total integration time with the 
SB3 bundle was 5400 s divided into three individual ex- 
posures of 1800 s each. The seeing conditions were fairly 
stable during the two periods of observation (about 0'.'7 
the first night, when we used the SB2 bundle, and I'.'l 
the second night when we used the SB3 bundle). 

The data reduction consists of three main steps: 1) the 
basic reduction of two dimensional fiber spectra, 2) the 
line fitting, and 3) the generation of maps of spectral fea- 
tures (such as intensity maps and velocity fields) from the 
calibrated spectra. The basic reduction was performed 



using the IRAF^ environment, and followed the stan- 
dard procedures applied to spectra obtained with two- 
dimensional fiber spectrographs (Arribas et al. 1997 and 
references therein). The absolute flux calibration was 
performed observing the standard star Feige 34 with the 
same instrument configuration and data reduction pro- 
cedures used on the galaxy. For each spectrum, the total 
flux, the radial velocities and the velocity dispersions (a) 
were measured by fitting single Gaussian functions to the 
observed emission line profiles using the DIP SO package 
(Howarth & Murray 1988). For the generation of the 
maps, a two dimensional interpolation was applied using 
the free software IDA (Garcfa-Lorenzo et al. 2002). The 
blue and red stellar continuum images were taken as an 
emission line- free filter centered at Arest/rame =4475 A 
and 6525 A respectively, and with a rectangular band- 
width of 150 A . The measures on the SB3 maps were 
done through an aperture of 3'.'0 in diameter, while for 
the SB2 maps the apertures diameter were I'.'O (IC 694) 
and 0'.'9 (NGC 3690, to avoid the overlap between re- 
gions). 

Three sample spectra, covering the entire used spec- 
tral range, are shown in Fig. |2| All the spectra show 
the strong (H/3, [Oiii]AA4959, 5007, Ha-H[Nii]AA6548, 
6584, HeiA6678, and [Sii]AA6716, 6731) and the weak 
(NIA5199, HeiA5876, [Oi]AA6300, 6364, and HciA6678) 
emission lines in different ratios, as well as interstellar 
absorption lines (NaiAA5890, 5896), that will be used to 
trace the neutral interstellar medium kinematics. 

Complementary HST multiwavelength images both 
in continuum and emission lines are also available. 
These additional data include the red (WFPC2-F814W, 
archive images) and the near-infrared (NICMOSl- 
FllOM, NICMOS2-F160W and NICMOS2-F222M, from 
AAHOO) continuum images. These filters are the HST 
equivalent to the ground-based Johnson-Cousins sys- 
tem /, J, H and K respectively (Origlia & Leitherer 
2000). Images of the ionized gas as traced by the Paa 
(Are^t=1.875 urn) light (NICMOS-F190N, from AAHOO) 
will be used too. The archival WFPC2 image was cali- 
brated on-the-fly with the best available reference files, 
while the NICMOS images were reduced with routines 
from the package NicRed (McLeod 1997; see details on 
AAHOO). 

3. SELECTION OF REGIONS UNDER STUDY 

In order to analyze the system we have selected sev- 
eral regions distributed over the individual galaxies, the 
interface between them and regions external to the galax- 
ies (Fig. nj. The selection includes compact unresolved 
high-surface brightness knots, star clusters and more dif- 
fuse, extended resolved low-surface brightness areas. For 
IC 694 and NGC 3690 we based our selection on that 
done by AAHOO, that includes the nuclear regions (A and 
Bl respectively) and the most important star- forming re- 
gions and bright star clusters. To identify the position 
of the regions under study in the lower resolution IN- 
TEGRAL maps, a direct comparison of the convolved 
WFPC2 image with the red continuum image was per- 
formed. The following sections present the main results 

* IRAF software is distributed by the National Optical Astron- 
omy Observatories (NOAO) which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) in coop- 
eration with the National Science Foundation. 
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obtained for these different regions and discuss their 
implications separately for the two individual galaxies, 
IC 694 and NGC 3690, and for the interface region be- 
tween them. The physical properties measured for each 
individual region are presented in Tables 1 to 3. 

4. STELLAR AND IONIZED GAS MORPHOLOGY 

Arp 299 is the result of a merger between the galax- 
ies IC 694 and NGC 3690, whose nuclei (A and Bl re- 
spectively, see Fig. ^ are located about 4.6 kpc apart. 
The system shows a morphologically disturbed structure 
characterized by the presence of compact high-surface 
knots, and a complex system of filaments and dust lanes, 
as already indicated by previous studies at several wave- 
lengths (e.g. Gerhz et al. 1983; Telesco et al. 1985; 
Baan et al. 1990; AAHOO). In what follows, we discuss 
the differences in the spatial distribution of the stellar 
and ionized gas components. 

4.1. Interface Region 

The two INTEGRAL/SB3 continuum images of 
Arp 299, obtained at wavelengths blueward of Ha and 
H/3, reproduce well the overall stellar morphology de- 
tected with the HST WFPC2 F814W image (Fig. EJl al- 
though at a lower spatial resolution (fiber diameter 2'.'7, 
scales of ~0.5 kpc). In spite of that, the interface region 
between the two galaxies, and the external regions show 
well defined features such as the bridge to the south of re- 
gion K5 connecting NGC 3690 and IC 694, and the elon- 
gated emission connecting Kl towards K12 source as part 
of a tail-like structure located to the SW of NGC 3690 
(see Fig. [TJ. 

The general structure of the low surface brightness ion- 
ized gas in the interface region, as traced by for example 
Ha (Fig. |3Jl, is similar to that of the stellar component. 
However, several differences exist. Diffuse ionized gas 
is detected NE of NGC 3690, associated with the K5 
and C knots; the latter region shows Paa emission and 
SN activity in radio (Neff et al. 2004), and is part of a 
bridge of ionized material (see Fig. |2| Ha map) to the 
north of the stellar bridge that connects the two individ- 
ual galaxies. The [Oiii]A5007 emission map presents an 
elongated structure located towards the south-east of Bl 
and B2 tracing an area where the ionized gas presents 
high excitation state (see H6.1|l . However, as showed in 
the next two sections, the most important morphological 
variations are associated to the main bodies of the galax- 
ies when observed with the higher angular resolution SB2 
bundle. 

4.2. IC 694 

Previous studies have concluded that IC 694 was a spi- 
ral galaxy that due to the merging process shows a more 
disrupted morphology (e.g., Gerhz et al. 1983; Standford 
& Wood 1989; Hibbard & Yun 1999). The stellar struc- 
ture as seen from both the HST and the INTEGRAL 
continuum images (see Fig. 0J) is divided into two main 
substructures: (1) The arm (region centered around A7) 
of about 1.8 kpc in length, that could be one of the origi- 
nal spiral arms disrupted and deformed as a consequence 
of the merging-process, and (2) the main body of the 
galaxy, formed by the dust-enshrouded infrared nucleus 
A (AAHOO), and surrounding regions. The two struc- 
tures are well separated (distances of ~1 kpc) by what 



seems to be a dust lane, that is clearly visible in the HST 
near-IR color maps (AAHOO). The effect of the internal 
extinction is remarkable in the continuum images as the 
relative surface brightness of the regions, in particular 
nucleus A and surrounding regions, changes with wave- 
length (see Fig. 0] and discussion in ii5.2|l . 

The ionized gas distribution (see Fig. shows sev- 
eral distinct features worth mentioning. The nucleus of 
the galaxy (A) appears as a weak optical line emitting 
source, whereas the regions of high surface brightness are 
located at different positions (e.g., part of the arm and 
towards east, and regions close to A) and separated up 
to 1 kpc from the nucleus. However, in the near-IR A is 
the brightest Paa an d Br 7 source (see AAHOO and Sugai 
et al. 1999, see also ^5.'2\i . indicating that the extinction 
towards the nucleus is very high (using Paa/Ha ratio, 
Ay=6.1±0.7). There are regions (e.g. knot G) that ap- 
pear as a faint continuum sources, but that are clearly 
identified as strong optical line emitters. In addition to 
differences in the morphology of the ionized gas (i.e. exci- 
tation conditions) with respect to that of the continuum, 
there are also differences in the ionized gas distribution 
as traced by several lines. The hydrogen recombination 
lines and the [Nii]A6584 and [Sli]AA6717, 6731 emission 
lines show a similar distribution, dominated by the emis- 
sion from regions located on the arm and towards the 
east, and in the regions surrounding the nucleus. How- 
ever, the shock tracer [Oi]A6300 line is concentrated in 
region HI, while the higher excitation line [Oiii]A5007 
appears to be brightest in regions A7 and G. These dif- 
ferences in the distribution of the emission line gas are 
not only due to internal extinction effects, but they also 
reflect changes in the ionization conditions as a function 
of location within the galaxy (see i l6.1|l . 

4.3. NGC 3690 

The stellar light distribution as observed with the HST 
shows an irregular morphology dominated by the emis- 
sion coming from the high surface brightness region B2 
(identified therefore as the optical nucleus), and from 
region Bl and complex C^ located at distances of 0.49 
and 1.3 kpc from B2, respectively (see Fig. EJ. However 
at near and mid-IR wavelengths the true nucleus of the 
galaxy (see discussion in i)10.2|l . Bl, becomes the bright- 
est continuum region in NGC 3690 (see AAHOO; Satyapal 
et al. 1999). This disordered morphology could be the 
combination of the original spiral structure with young 
star- forming regions such as the complex C, and regions 
D3 or C, newly formed as a consequence of the interac- 
tion process as suggested by several authors (e.g. Gerhz 
et al. 1983; Standford & Wood 1989; Baan & Haschick 
1990; Hibbard & Yun 1999). 

The distribution of the ionized gas (Fig. O is sim- 
pler than that of IC 694. The optical nucleus (B2) ap- 
pears as an extremely faint line emitter, and the struc- 
ture of the ionized gas is dominated by two high sur- 
face brightness regions separated by 1.8 kpc, nucleus Bl 
and the complex C (with the peak located towards C2), 
as previously founded with near-IR hydrogen recombi- 
nation lines (Sugai et al. 1999; AAHOO). The bright- 

^ As explained in Note 7, we will use the term complex C intro- 
duced by Gehrz et al. (1983) for the regions CI, 02 and C3 (see 
Fig. 0. 
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est emission line peak is associated to C2 for all lines 
but the shock tracer [Ol]A6300. For the latter, the nu- 
cleus Bl is the brightest emitter, as also observed with 
the near-IR [Fe ii]A1.64 /xm line (AAHOO). Analyzing 
in detail the secondary emission line peaks associated to 
Bl and surrounding areas, it is interesting to note that 
the high excitation [O iii]A5007 emission is dominated by 
Bl, while the Ha and H/3 emission peaks are associated 
to region B16 (which is located 0.26 kpc southwest of 
Bl, and identified in the A4475 continuum map). Low 
excitation coUisional lines [Nii]A6584, [Sii]AA6717, 6731 
appear to have secondary peaks in both Bl and B16. 
These changes in the emission line distribution repre- 
sent changes in the ionization conditions, and most likely 
in the ionizing source (AGN, young stellar clusters and 
shocks) on scales of less than 0.3 kpc. 

5. INTERNAL DUST DISTRIBUTION AND EXTINCTION 
EFFECTS 

As mentioned in the previous section, part of the com- 
plex structure observed in the continuum and emission 
line light distributions is due to large amounts of dust 
producing spatially dependent internal extinction effects. 
The two dimensional extinction maps of the gas (Figs. IHI 
[3 and (SI have been derived using the Ha/H/3 line ratio, 
a foreground dust screen model, and a mean interstellar 
extinction law based on that of Savage & Mathis (1979). 
No correction for the presence of underlying stellar hy- 
drogen absorption lines has been applied. Assuming an 
equivalent width value of EWabs (Ha) = EWabs (H/3) ~ 
2 A , the derived extinction could be overestimated by 
as much as AE(B-V) ^ 0.6 for regions with extremely 
faint H/3 emission (e.g. EW(H/3) < 20 A). The values 
for the visual extinction (Ay=3.1xE(B-V)) in particu- 
lar regions are given in Tables ^ |21 and 13 The results 
suggest that the highest extinguished regions are, in gen- 
eral, associated to the nuclei and some particular regions 
(like C), while the external Hii regions are mostly less 
affected by the presence of dust. 

5.1. Interface Region 

The extinction in the whole area covered by the SB3 
bundle on scales of about 3'.'0 (i.e. 0.6 kpc), is relatively 
uniform (see Fig. EJ. It has an average value of Ay = 1.9 
mag on the selected individual regions, and ranges from 
Ay(K2) = 0.9±0.7 mag to Ay(K5) = 2.9±0.7 mag (see 
Tabled for specific values on all the individual regions). 

5.2. IC 694 

On scales of about I'.'O (i.e. 0.21 kpc), the internal ex- 
tinction in the visual covers a wide range of values from 
about Av=0.6±0.5 mag in the southwestern section of 
the galaxy to Ay=3.4±0.5 mag in region H2, close to 
the nucleus A (see Table|21for specific values, and Fig. [3 
for the two dimensional E(B-V) distribution). As can be 
seen from Fig. |7| the Ha extinction corrected light dis- 
tribution shows an overall structure similar to that of the 
Paa line except for the region associated with the peak 
of the emission. While the nucleus A, identified as the 
brightest near-infrared source in the HST images (Fig. 
see also AAHOO), is also the dominant Paa line emitter, 
the peak of the extinction corrected Ha distribution is 
located in regions HI and H2, 0.3 and 0.5 kpc towards 



the east of the nucleus. Also region F, located 1.58 kpc 
to the south-east of A, appears as a strongly absorbed 
region (Av=3.2±0.5 mag), and as a strong Ha emitter 
after correcting for extinction. Additional secondary Ha 
sources are mostly associated with star-forming regions 
(the arm, a region towards the west of A) at distances of 
1 to 1.6 kpc from the nucleus, and delineating what seems 
to be a spiral-like arm structure previously identified in 
the Paa emission line map. 

5.3. NGC 3690 

This galaxy presents the widest range of extinction (de- 
rived from the Ha/H/3 ratio) for the entire Arp 299 sys- 
tem with values that for the individual regions range from 
Ay(C4)=0.5±0.5 mag to Ay (B2)=3.9±0.6 mag (see Ta- 
ble OJ. Also, the nucleus Bl in this galaxy is the region 
of the Arp 299 system which has associated the highest 
Ha/H/3-based extinction. Ay ~4.6 mag (Fig. |HJ). 

The overall structure of the extinction corrected Ha 
line distribution agrees well with that of the Paa line 
(Fig. |S1 see also AAHOO) not only on the high surface 
brightness regions such as the nucleus Bl and CI, but 
also on the low surface brightness region, including B2 
identified as the optical nucleus. 

5.4. Extinction estimate using optical and near-IR 
emission lines 

Additional line ratios involving IR hydrogen recombi- 
nation lines (e.g., Paa/Ha, Br7/Paa, and Paa/Pa/3), 
can provide better estimates of the internal extinction in 
highly extinguished regions as those identified in Arp 299 
(see Table EJ. 

The Paa image of the individual galaxies (see Figs. 
IHl and AAHOO) has been used for measuring all the 
regions under study with apertures of I'.'O for IC 694 
and of G'.'9 for NGC 3690 to match our ground-based 
apertures; the INTEGRAL Ha fluxes are given in Ta- 
bles El and 13 In most regions extinction estimates based 
on only optical and optical+infrared lines agree within 
the uncertainties (Table 0J. On the other hand, there 
are two regions, the nucleus of IC 694 (A) and C, 
where the use of the optical-|-infrared lines indicates a 
larger extinction than using only the optical lines, i.e. 
Ay (Paa/Ha) ~2x Ay (Ha/H/3). This means that the ex- 
tinctions derived using the Ha/H/3 ratio are in general 
consistent with the ones derived using the Paa/Ha ratio, 
except in highly absorbed regions, where there is an un- 
derestimate when only the optical lines are used. These 
extinctions are similar with those measured in the M51 
nuclear Hii regions (Calzetti et al. 2005), and in other 
LIRGs (Alonso-Herrero et al. 2006) using the same near- 
IR/optical emission lines. 

6. TWO-DIMENSIONAL IONIZATION STRUCTURE 

The ionization level for the individual regions identified 
in Arp 299 has been derived using standard optical emis- 
sion line diagnostic diagrams (BPT diagrams: Baldwin, 
Phillips & Terlevich 1981; Veilleux & Osterbrock 1987). 
H ii-like ionization (i.e. star formation) is present on the 
majority of the regions of the individual galaxies, except 
for nucleus Bl which clearly has Seyfert-like (i.e. pres- 
ence of an AGN) state. LINER-likc ionization is mainly 
present on the interface (Fig. |51 see also Tables 1,2, and 
3 for specific values of the emission line ratios). Several 
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regions slrow clranges in the classification depending on 
the specific diagram used; these variations have been in- 
terpreted as a consequence of changes in the dominant 
ionizing mechanism from photoionization to shocks (e.g. 
Monreal-Ibero et al 2006 and references therein). In par- 
ticular, the [Oiii]A5007/H/3-[Oi]A6300/Ha line ratio is 
more sensitive, and therefore a more reliable tracer, to 
changes in the ionization conditions due to fast shocks 
than any of the other optical emission line ratios in- 
volving [Nii]A6584/Ha or [S ii]A6725/Ha ratios (Dopita 
& Sutherland 1995). In the following subsections a 
more detailed discussion of the two-dimensional ioniza- 
tion structure of the ionized gas distribution based on 
the [Oiii]A5007/H/3-[Oi]A6300/Ha diagram and on the 
[Oiii]A5007/H/3, and [Oi]A6300/Ha maps for the inter- 
face, IC 694, and NGC 3690 is given. 

6.1. Interface region: Evidence for an off-nuclear 
Seyfert-like ionization cone 

The [Oiii]A5007/H/3-[Oi]A6300/Ha line ratios have 
been obtained for aU the INTEGRAL/SB3 individual 
spectra (Fig. [TUIl . Although these data include the 
high surface brightness nuclear regions of IC 694 and 
NGC 3690, and their external regions, they also cover 
the area that has been identified as the interface zone, 
i.e. the extended low surface brightness region connect- 
ing the two galaxies up to distances of about 4 kpc from 
their nuclei. Changes in the ionization structure traced 
by the [O iil]A5007/H^ and [Oi]A6300/Ha two dimen- 
sional spatial distribution line ratios are also shown (Fig. 

cni). 

Three different types of ionization, Hii-, LINER- and 
Seyfert-like, are identified, and the line ratio maps clearly 
show that the ionization structure has a well defined spa- 
tial distribution. The majority of the regions associated 
with the two individual galaxies (marked with two iso- 
contours in the figure) are dominated by H ii-like ioniza- 
tion, as expected from the intense star formation taking 
place there (see more details in Sj^. There are some 
exceptions; the nucleus Bl shows LINER- like activity 
at this angular resolution. The Seyfert-like ionization 
is clearly resolved with a size of about 7" (^1.5 kpc) 
and has a conical morphology, with an opening angle 
of about 54° and its apex located in nucleus Bl (at a 
projected distance of about 1.5 kpc from the peak of 
the nebula) within the angular resolution (see Fig llU|) . 
Moreover, this highly ionized gas is not associated to 
any particular concentration of stellar mass according to 
the optical (WFPC2/F814W; see Fig. [TJ HST contin- 
uum image. Off-nuclear Seyfert-like nebulae at distances 
of few kpc from the nucleus have already been reported 
in other (U)LIRGs with Seyfert nucleus (Mrk 273: Col- 
ina et al. 1999), and interpreted as photoionization of 
extranuclear interstellar gas by the AGN. Therefore, our 
data strongly suggest photoionization by a radiation cone 
escaping from a central dust-enshrouded AGN source 
located in Bl. The highly ionized conical structure in 
Arp 299 (NGC 3690) is also detected at distances of up 
to 4 kpc from Bl; however, the ionization associated to 
these outer regions of the cone is LINER-like rather than 
Seyfert (Fig. 110(1 . Therefore, it seems that within the 
cone there are two ionization regimes well separated spa- 
tially. At projected distances from Bl smaller than about 
2 kpc the ionization is likely due to radiation coming from 



the AGN located in Bl, while at larger distances LINER- 
like ionization is dominant. This LINER-like ionization 
is also found in other areas of the interface region, how- 
ever with a lower degree of excitation, as traced by the 
[O i]A6300/IIq! ratio measured in these regions. 

6.2. IC 694 

Over the area covered by the INTEGRAL SB2 bun- 
dle, Hii and LINER-like ionization are present in this 
galaxy (Fig. 111(1 . The main body of the galaxy, includ- 
ing the regions close to the nucleus A and the spiral 
arm structure (A7, G) towards the east (F), is dom- 
inated by the star- forming knots. LINER-like ioniza- 
tion is present elsewhere, being more prominent in dif- 
fuse low surface brightness regions located at a projected 
distance of about 1.8 kpc. The [N ii]A6584/IIa and the 
[S ii]A6725/Ha line ratio maps (not shown) also have a 
similar structure though, as explained above, they tend 
to indicate a lower ionization. 

6.3. NGC 3690 

Three types of excitation level (Hii-, LINER- and 
Seyfert-like) are present in this individual galaxy, as 
shown in the [O iii]A5007/H/3-[O i]A6300/Ha diagnostic 
diagram and in the two dimensional maps (Fig. 112(1 . 
Regarding the nucleus Bl, these higher resolution data 
give a clear Seyfert-like classification in all these line ra- 
tios (Fig. El, and show that there is an extended zone 
of about 2" (0.4 kpc) in length associated to Bl with 
Seyfert characteristic. This is the inner section of the 
AGN-like ionization cone detected with the SB3 bundle, 
with its apex coincident with Bl. Moreover, the region 
B2 presents also a high ionization state in the border line 
of Seyfert/LINER, characteristic that could indicate the 
presence of the counter ionization cone. The detection of 
Seyfert-like activity based on optical emission lines rep- 
resent the first evidence in the optical for the presence 
of an AGN in Bl, and supports the recent claim of an 
AGN based on X-ray (Delia Ceca et al. 2002) and mid-IR 
(Gallais et al. 2004) data. Previous spectroscopic stud- 
ies in the optical (Gehrz et al. 1983; Coziol et al. 1998; 
Keel 1984) did not find AGN emission, probably due to 
the large size of the slits (3" 6, 2'.' 5, and 8'.'0). This is 
in agreement with our results, where the Seyfert nature 
of the Bl nucleus has been diluted in the SB3 {0— 2'! 7) 
spectra, whereas in the higher spatial resolution SB2 {0— 
0'.'9) spectra the classification is clear. 

H ii-like ionization dominates the rest of the main body 
of the galaxy (e.g. complex C), while LINER-like is 
present in a few outer regions. The [O i]A6300/Hq! map 
has its peak displaced about 0'.'7 (0.14 kpc) to the north 
with respect to Bl. Similar displacements are observed 
in the [N ii]A6584/IIa and in the [S ii]A6725/IIa maps. 
These lines and in particular [O i]A6300/IIa are a good 
tracers of shocks, thus these shifts may indicate differ- 
ent ionization mechanisms around Bl on scales of 200 
pc. These small scale ionization changes need further 
investigation with high angular resolution spectroscopy. 

7. ELECTRON DENSITY 

The two-dimensional electron density maps of the ion- 
ized gas distribution have been derived for the entire 
Arp 299 system using the standard [S ii]A6716/[S ii]A6731 
lines ratio (AUer 1984; Osterbrock 1989). The average 
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line ratio for the interface and the individual galaxies 
show a similar value, about 1.3±0.1, indicating an aver- 
age uniform electron density (~100 cm~'^) over the entire 
Arp 299 system. However, the ionization cone detected 
in NGC 3690 presents a well defined gradient of this line 
ratio value, that roughly ranges from 1.2±0.15 in the 
higher excitation (i.e. Seyfert-like) center of the cone, to 
1.8±0.15 in the lower excitation (i.e. LINER-like) outer 
section. Thus, the electron density seems to decrease 
from about 180 cm~^ at a distance of 1.5 kpc from Bl 
to less than 10 cm^'^ at about 4 kpc. 

8. THE ORIGIN OF THE LINER-LIKE IONIZATION 

Most regions of the Arp 299 system have been clas- 
sified as LINER-like, although the mechanism produc- 
ing this ionization needs to be determined. Some au- 
thors have found in (U)LIRGs a positive correlation be- 
tween the [O i]A6300/Ha ratio and the velocity disper- 
sion (Monreal-Ibero et al. 2006), which has been con- 
sidered as evidence supporting the shocks as the origin 
for the ionization (Armus et al. 1989; Dopita & Suther- 
land 1995). The [Oi]A6300/Ha versus the Ha velocity 
dispersion for all the spectra obtained with INTEGRAL 
has been represented (Fig. I13|) in order to investigate 
the origin of the LINER-like ionization in Arp 299. The 
four main regions of the system observed with the IN- 
TEGRAL SB2 and SB3 bundles are analyzed separately: 
the inner and outer sections of the ionization cone, the 
interface region and finally the individual galaxies. The 
outer section of the ionization cone (stars in Fig. ^1 
Left panel. LINER-like excitation , i.e.[Oi]A6300/Ha 
> -1.25) is independent of the velocity dispersion and 
present a stable ionization level, while with a wider va- 
riety in the ionization level, the inner section of the 
cone (Seyfert-like, triangles in Fig. ^1 left panel) does 
not follow the correlation founded in external regions of 
ULIRGs (Monreal-Ibero et al. 2006). The interface re- 
gion (crosses in Fig. E| Left panel) present a stable 
LINER-like excitation level, and is not dependent of the 
velocity dispersion. 

With regards to the individual galaxies IC 694 and 
NGC 3690 as observed with the lower angular resolution 
(SB3) bundle (filled circles in Fig. ^1 Left panel), the 
ionization is mainly Hll-like, with no evidence of a kine- 
matical correlation. For the higher angular resolution 
(SB2) data (see Fig. ^1 Center and Right panels), the 
actual measurements of the individual galaxies do not 
follow the correlation, although they present a well de- 
fined clustering. The galaxy IC 694 has a mean log(T~1.8, 
and presents about the same fraction of H ii and LINER 
activity. In the case of NGC 3690 the velocity dispersion 
is clustered around logcr^2, and there is a large fraction 
of regions with H ii-like ionization. 

Clearly the present data show a different behavior than 
that found in ULIRGs (Monreal-Ibero et al. 2006), and 
this is probably due to the fact that on each case the 
measurements are done over regions located at different 
distances to the nucleus, and with different character- 
istics. The original correlation was derived for diffuse 
extra-nuclear regions in ULIRGs, which are located at 
distances of about 5—15 kpc from the nucleus of the 
galaxy. In the present case the interface region and the 
ionization cone have regions at distances of up to 4 kpc 
from the nucleus, and comparing with the original lim- 



its of the correlation (see shaded area on Fig. ^1 Left 
panel) the data of the interface and both regions of the 
ionization cone lie in general above the ratio, i.e. this 
regions show higher excitation for the same velocity dis- 
persion. In the case of the individual galaxies observed 
with the actual higher resolution bundle (Fig. 1131 Center 
and Right panels) the present data are closer to the nu- 
cleus than in the case of ULIRGs (up to '^2.5 kpc) and 
in most cases the regions involved are compact bright 
sources rather than diffuse. Comparing with the original 
limits of the relation (see shaded area on Fig. Center 
and Right panels) the IC 694 data are clustered around, 
while the ones of NGC 3690 are below the correlation. 
This means that in the present case for the individual 
galaxies the star formation influence is still present, and 
that the velocity dispersion close to the nucleus traces 
mass rather than flows/shocks. 

Given an ionization mechanism, the ionization state 
of the surrounding gas can be changed as a result of 
the variations in the electron density, that affect to the 
ionization parameter ( [/ocNe~^r~^). This appears to be 
the case in the outer high excitation cone associated to 
NGC 3690. As already mentioned, the electron density 
seems to have a radial variation within the cone, and 
by having the exact variation of Ng with r, the radial 
dependence of the parameter U would be described more 
precisely. Unfortunately in the present case the exact 
radial dependence of Ne cannot be derived; in spite of 
that, the variation detected on the electron density, and 
therefore the LINER activity could be accounted for by 
the influence of an AGN (if Ng oc r~"). The variation of 
Ng implies a variation on the ionization parameter (U (x 
r~^+"), and thus the low excitation is valid if the source 
is not local. 

In conclusion, for the Arp 299 system studied here 
there are no clear evidences for tidally induced shocks as 
found in ULIRGs. Therefore, the origin for the LINER 
ionization is unclear, and at least in some zones (i.e. 
outer parts of the cone) it could, perhaps, be due to a 
radiation field similar to that of a low-luminosity AGN. 

9. STAR FORMATION PROPERTIES 

It is well known that merging processes between gas 
rich galaxies can trigger star formation. As predicted 
by numerical simulations (Mihos & Hernquist 1994), the 
process can vary depending on whether or not the parent 
disk galaxies have a bulge. The star-forming process at a 
low level (SFR~20 Mq yr~^) consequence of the merging 
of two disk galaxies without a bulge is produced early on 
in the merging phase, and can be sustained for about 
150 Myr. On the other hand, in the merger of galaxies 
with a bulge the starburst (SFR ~ 50-100 Mq yr'^) is 
produced later in the merging phase, and lasts for about 
50 Myr. 

To investigate the properties of the young stellar pop- 
ulations of a galaxy assuming instantaneous and contin- 
uous bursts, synthetic models such as SB99 (Leitherer et 
al. 1999; Vazquez & Leitherer 2005) can be used. In the 
following the stellar population synthesis model consid- 
ers a Salpeter initial mass function (IMF), with lower and 
upper mass limits of 0.1 M0 and 100 Mq, the Geneva 
tracks with high mass loss, and the Calzetti extinction 
law (Calzetti et al. 2000). Instantaneous bursts are nor- 
malized to 10^ Mq, while the continuous star formation is 
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normalized to a star formation rate of 1 Mq yr~^. Com- 
bining the information from the two dimensional distri- 
bution of the Ha equivalent width and luminosity, and 
the photometric results, parameters such as the extinc- 
tion to the gas, the age of the burst, and the stellar mass 
for each individual region can be derived. For the present 
analysis we use the INTEGRAL SB2 bundle with a reso- 
lution of 0'.'9 (i.e. 0.18 kpc), thus in the majority of cases 
the data include both the contribution from the stellar 
clusters and from a diffuse component. 

9.1. Photometry of the regions of IC 694 and NGC 3690 

The aperture photometry as derived by AAHOO, us- 
ing the HST NICl FllOM, NIC2 F160W, NIC2 F222M 
images is used in the following analysis. In addition opti- 
cal (WFPC2/F814W) photometry over the same regions, 
using the same apertures, and after performing the back- 
ground subtraction and the aperture correction (for the 
latter see Holtzman et al. 1995) has been obtained. The 
apertures used were 0'.'9 diameter for the nuclei of IC 694 
and NGC 3690 A and Bl, and for the star-forming re- 
gion B2, and 0'.'5 diameter for the rest of the star- forming 
regions. The HST magnitudes and colors derived are 
shown in Table [S] In addition, the HST magnitudes and 
colors were modeled by convolving the simulated galaxy 
spectral energy distribution with the system (optics, fil- 
ter and detector) transmission hmction. 

9.2. Properties of the Ha emitting regions 

In the following, instantaneous bursts models charac- 
terized by the predicted Ha luminosity and equivalent 
width, and near-IR colors will be used to study the stel- 
lar populations of the Ha emitting regions under anal- 
ysis. Comparing the measured values with the models, 
it is possible to derive the age, extinction and mass of 
the individual regions under study. Regarding to the nu- 
cleus Bl, and as previously indicated, it is classified as 
Seyfert, and the fraction of the ionizing flux coming from 
the buried AGN in this region is not well determined. 
Therefore the conclusions regarding the star formation 
properties in Bl based on the Ha luminosity have to 
be taken with caution, and represent an upper limit. Be- 
cause of this fact, the nucleus Bl is basically not included 
in the following discussion, although the properties de- 
rived are listed in the tables. 

The Ha luminosity in specific regions with large Ha 
equivalent width (Fig. I14() traces preferentially the pres- 
ence of young stellar populations (<5— 10 Myr). How- 
ever, the presence of an old underlying population tends 
to decrease the EW, and therefore to underestimate the 
age of the young stellar population. In the case of IC 694 
and NGC 3690 the respective nuclei A and Bl, and 
the NGC 3690 brightest optical region B2 present low 
EW(Ha). The ages derived from the EW(Ha) for sev- 
eral individual regions are shown in Table El 

With an age for the nucleus (A) of 5.3x10® years, 
IC 694 does not show a wide spread in age for the indi- 
vidual star-forming regions; almost all are around 5x10® 
years, except for II with 5.9x10® years and the youngest 
regions, concentrated towards the east of the spiral-like 
arm (F and G, with 3.3 and 3.4x10® years respectively). 
The galaxy NGC 3690 shows an average age for the in- 
dividual regions of ^5.5x10® years, with the eldest one 
(7.2x10® years) located towards the north of B2. The 



youngest sources are located in the northern part of the 
galaxy, and correspond to complex C. Independent esti- 
mates (Satyapal et al. 1999) also based on instantaneous 
bursts and characterized by a Salpeter IMF, derived ages 
of about 6x10® years for the main components of the 
Arp 299 system, that are compatible with the results 
obtained here. Compared to AAHOO, who used a dif- 
ferent model with Gaussian bursts, the present age for 
B2 is younger, but applying their EW(Br7)<5 A to our 
model, the age for B2 is >7xl0® years, compatible with 
the actual calculations. 

Extinction and dust properties of the star-forming re- 
gions can be studied using the near-IR colors FllOM- 
F160W (approximately J-H) and F160W-F222M (ap- 
proximately H-K) together with the EW(Ha). Under 
the assumption that the extinction to the stars is equiva- 
lent to the extinction towards the gas, the EW(Ha) value 
is independent of Ay fFigll5(l . We have unreddened both 
observed colors using a foreground screen model and the 
extinction derived from the Ha/H/3 ratio, and compared 
with the models. Once corrected for extinction, and tak- 
ing into account the error bars, in the majority of the re- 
gions the values are close to the predicted Ay=2— 3 mag 
track. Thus in general the extinction derived from the 
Ha/H/3 ratio presented in Tables |2] and 13 show an offset 
with respect to the ones predicted by the colors. There 
are two possible explanations for this discrepancy; on the 
one hand the extinction derived from the Ha/H/3 ratio 
can be underestimated (see ii5.4ll . on the other hand the 
discrepancy traces the contribution of hot dust emission. 
For the highly reddened regions (e.g. A, Bl^*', CI) the 
colors measured (and corrected by extinction) are very 
red. As already mentioned (and discussed by AAHOO) 
one of the most plausible explanations for this color ex- 
cess is the presence of hot dust (T>600 K), which would 
mainly affect the F160W-F222M color. In an attempt 
to take into account the contribution of the hot dust 
emission, the A~^-Ba wavelength modified black-body 
thermal emission (Hildebrand 1983) has been used. The 
results demonstrate that depending on the temperature 
(which has not necessary to be homogeneous over the 
line of sight) and the amount of dust, the contribution of 
hot dust (~800 K) emission could represent, on average, 
about 80% of the K-band emission (Fig. El Right panel) 
for the highly extincted regions (A, Bl, CI), and about 
30% for the rest of the regions. This result supports our 
previous conclusion that the dust is more concentrated 
towards the nuclei and some particular regions, than in 
the extra-nuclear star forming regions. 

An alternative to the presence of a hot dust emission 
is that the structure of the dust-gas distribution does 
not follow the simple screen geometry, but it is closer 
to a clumpy foreground or a mixed dust-star distribu- 
tion. Assuming a model where the gas is homogeneously 
mixed with the stars, the effect of an Ay=6 mag over the 
F160W-F222M and F110M-F160W measured colors wiU 
produce bluer colors: about 0.3 and 0.4 mag respectively 
(Scoville et al. 2000; see Fig. I15|l . For the extreme case 
Ay=30, the corresponding shifts will be 0.6 and 0.8 mag, 
i.e. in region CI the F160W-F222M color will move from 

In this discussion related to the hot dust the nucleus Bl is 
included, because even with high errors in the age derived from the 
EW(Ho), the color excess exists and needs to be explained. 
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the observed 2.67 to 2.07. Thus the colors modeled with 
this alternative dust geometry would imply that a hot 
dust contribution is not mandatory in the external less 
extincted regions, but it is still required to explain the 
color excess, in the highly extincted regions (A, Bl and 
CI). 

The stellar mass of the young population has been 
measured using two different indicators, the Ha luminos- 
ity and the MigpvF absolute magnitude, and comparing 
with the predicted values for the ages derived using the 
EW(Ha). In both cases the Ha/H/3 ratio has been used 
for correct by extinction (for A and C the Paa/Ha ratio 
was used too). 

The results obtained using Migoiv (see Fig. ^1 and 
Table show that the regions with more young stellar 
mass content are nucleus A and regions B2, B21 and 
CI (460, 300, 58 and 81x10^ Mq respectively). The 
stellar mass content for the rest of the regions, including 
the low-surface brightness knots, ranges between 3.1 and 
14.0xlO^Mo (details on Table |SJ). Previous derivations 
of the stellar mass of the clusters using absolute H-band 
magnitudes vs H-K (not corrected for extinction) with a 
different stellar populations synthesis model range, from 
about I.SxIO^Mq to S.OxlO^Mg (AAHOO). 

Using the Ha luminosity, the young stellar mass con- 
tent is again concentrated on the nucleus A, and regions 
B2, B21 and CI (655, 679, 74 and 74x10^ Mq respec- 
tively). For region C, where there is not near-IR pho- 
tometry, the mass derived using the Ha luminosity is 
115xlO^Mo. The rest of the derived masses range from 
5.7 to 57x 10^ Mq. Comparing the mass in newly formed 
stars with the derivations of AAHOO, there are no signifi- 
cant differences except for B2 where we derive a value two 
times higher. The masses of the young Ha-emitting re- 
gions identified in Arp 299 are within the range of masses 
founded in LIRGs (3x1O6-1O8M0, ScoviUe et al. 2000; 
Alonso-Herrero et al. 2002) and are on average about 
100 times more massive than the galactic globular clus- 
ters (IO^Mq, van den Bergh 1995), probably due to the 
fact that the present regions have conglomerates of less 
massive clusters. 

When comparing the stellar masses derived with the 
Ha luminosity with the mass obtained from the abso- 
lute F160W magnitude, there is a good agreement in the 
nucleus A, but for the rest of the regions there are differ- 
ences from a factor '^l to 5. One of the main reasons for 
this difference is that in the majority of the cases (except 
A and B2) an aperture of 0'.'5 has been used for the HST 
cluster photometry, instead the 0'.'9 or I'.'O used for the 
INTEGRAL data analysis (for the EW(Ha) and L(Ha) 
estimates); this difference can lead to offsets of a factor 
3 in the mass estimates. Factoring in the aperture ef- 
fects, the two methods give the same mean value within 
a factor 2. 

9.3. Contribution to the bolometric luminosity 

The derived ages of the young Ha-emitting star- 
forming regions cover a range of between 3x10^ years 
to at least 7-8x10^ years old. This is consistent with 
the methodology used here: since the Ha traces prefer- 
entially the massive young ionizing stellar clusters, pop- 
ulation older than 10 Myr is difhcult to be detected as 
Ha decreases. In any case there is further evidence that 
the star formation in Arp 299 has persisted for at least 



the last 10-15x10'' years (AAHOO), with high activity on 
the star-forming regions. On the basis of this dispersion, 
we consider that the overall star formation of Arp 299 
can be modeled by a continuous star formation. Un- 
der this premise the star formation rate (SFR) and the 
IR luminosity contribution can be derived from the ex- 
tinction corrected Ha luminosity, using Kennicutt (1998) 
expressions, valid for a Salpeter IMF with mass cutoffs 
of 0.1-100 Mq: 

SFR{Mq year^^) = 7.9 x lQ-'^'^L{Ha){erg s'^) (1) 

Lmierg s'^) = 175.5 x L{Ha){erg s'^) (2) 

that strictly speaking are only valid for starbursts with 
ages less than 10^ years. 

Without taking into account the nucleus Bl (AGN can- 
didate), the SFR of the entire system is about 43 M0 
year"-'^, and the Lir is 9.5x10"*^ erg s~^, whereas the 
luminosity derived from the IRAS fluxes is 2.2x10^^ erg 
s"^, thus the contribution of the studied regions to the 
IR luminosity is about 45%. There is a lack of about 55% 
in luminosity, that can be accounted for by the AGN lo- 
cated in Bl. Previous estimates based on the mid and 
far-IR spectra of Arp 299 assign a 20% of the Lm com- 
ing from Bl and B2 (Charmandaris et al. 2002), while at 
11.7 /im the percentage of flux coming from Bl and B2 is 
about 30% (Keto et al. 1997). So even considering that 
the majority of this contribution comes from Bl, there is 
still a lack of about 25—35%. This inconsistence is prob- 
ably due to the fact that the present result is limited by 
the derived extinction and/or the models. 

10. 2D GAS KINEMATICS AND DYNAMICAL MASS 
DERIVATIONS 

Previous studies of Arp 299 (Gehrz et al. 1983; Tele- 
sco et al. 1985) concluded that this system consists 
of two spiral galaxies disturbed by the merging pro- 
cess, which are rotating in opposite sense (Augarde & 
Lequeux 1985). More recently Hibbard & Yun (1999) 
have described Arp 299 as a structure formed after a 
prograde-retrograde encounter between IC 694 (Sab-Sb) 
and NGC 3690 (Sbc-Sc) that took place some 750 Myr 
ago, and having a long tidal tail (180 kpc in length) de- 
tected in Hi. Although there are still two structures 
clearly distinguishable with the two nuclei separated by 
about 4.6 kpc, the warm gas component has lost its kine- 
matic identity and the velocity field presents a very ir- 
regular pattern (Rampazzo et al. 2005). 

10.1. Gas kinematics 

The kinematics of the ionized and the neutral gas is 
traced by the Ha emission line and the NaiA5890 ab- 
sorption line respectively (Sparks et al. 1997). The ion- 
ized gas mean velocities and velocity dispersions were 
obtained by fitting the Ha+[N ii]AA6548,6584 complex 
with a single Gaussian velocity component per line; the 
same method was applied to the He iA5876-l-NaiAA5890, 
5896 complex. The values derived for the individual re- 
gions are presented in Tables |H1 and El For the veloc- 
ity measurements, the errors have been calculated tak- 
ing into account both the systematic uncertainties of the 
wavelength calibration (same for all lines) in the velocity 



Integral field spectroscopy of Arp 299 



9 



measureraent, and tlie uncertainties in tlie measurement 
of the centroid of the line (that depends on the S/N of 
the studied line). One of the most relevant results ob- 
tained is that the two dimensional structure of the ion- 
ized and neutral gas velocity field show similar features 
in both the individual galaxies and the interface region, 
indicating that the ionized and neutral gas kinematics 
are coupled. 

The structure of the Ha velocity dispersion shows that 
there is a variation on the average velocity dispersion 
values: the interface region has the lowest (45 km s^^), 
IC 694 the intermediate (65 km s'^) and NGC 3690 (120 
km s~^) the maximum values of a. 

10.1.1. Interface Region 

The overall velocity field of the ionized gas as observed 
with the SB3 bundle has a complex and irregular pat- 
tern, that does not indicate the presence of virialized 
simple (e.g., rotating) motions, on scales of 0.6 kpc or 
larger (Fig. El Left). The structure of this map agrees 
well with that obtained by Rampazzo et al. (2005) cov- 
ering a larger area (FoV of 5'. 9) on the sky. The veloci- 
ties in the northern part are receding, and in the south- 
ern part are approaching to the systemic value, defined 
as the velocity in nucleus A (the dust-enshrouded nu- 
cleus of IC 694) with the higher spatial resolution bun- 
dle SB2, Va=3121±27 km s'^ (see section i JlO.l.^^ . The 
corresponding lower spatial resolution SB3 velocity is 
Va=3190±27 km s^^. 

The velocity field of the neutral gas obtained from the 
NaD absorption lines, despite of being noisier (the in- 
tensity of the NaD can be typically 2% of that of Ha), 
shares similar structural characteristics with the ionized 
gas velocity field (Figure El Center panel), and presents 
a systemic velocity of 3057±34 km s~^. 

The ionized gas velocity field presents a peak-to-peak 
velocity of about 230 km s~^, while for the neutral gas 
the variation is smaller, about 170 km s""'^. In both cases 
the area with the largest approaching velocities is close to 
B2. The ionization cone does not show relevant features 
in the ionized gas velocity field, nor in the neutral gas 
velocity field. With the exception of the nuclear regions, 
the ionized gas in the galaxies does not have a distinct 
kinematics with respect to the interface region (i.e., the 
structure of the velocity field does not correlate with the 
distribution of mass as traced by the continuum images). 

The velocity gradient of the ionized gas through the 
interface is about 40 km s~^ kpc~^, compatible with 
the Ha velocity gradients previously measured across 
ULIRGs observed with a slit (up to 47 km s~^ kpc~^ 
for galaxies with double nuclei, Martin 2006). 

The velocity dispersion shows distinct features over the 
entire area; there are large differences in cr, that allow us 
to discriminate spatially between quiet and turbulent gas 
with macroscopic motions (Fig. ^\Right panel). The so- 
called interface has an average cr value of 45 km that 
traces no turbulent gas. The low velocity amplitude and 
velocity dispersion of the interface indicates that the ion- 
ized gas could be in a close to quiescent phase or slowly 
rotating. Larger a values are measured on the previously 
defined ionization cone, with the maximum (160±21 km 
s^^) located on the peak of Seyfert-likc ionization. This 
spatial coincidence could indicate that in that area the 
two contributions of the biconic surface of the cone are 



added, and thus the line measured is wider (Heckman 
et al. 1990). On the other hand, the maximum veloc- 
ity dispersion would imply that the largest macroscopic 
motions, with approaching velocities, are taking place in 
that highly ionized area. 

10.1.2. IC 694 

The velocity field of both the ionized and neutral gas 
presents a similar structure (Fig. 1181 Left and Center 
panels), when observed with the SB2 bundle. It has 
an irregular structure, with the largest receding veloc- 
ities associated with region I (3181±27 km s~^ for the 
ionized gas and 3155±34 km s~^ for the neutral), lo- 
cated at 0.42 kpc SW from the nucleus (A). These maps 
show a general trend with velocities changing from NW 
(receding) to SE (approaching), in agreement with the 
pattern found with the SB3 data in this region. The 
systemic velocity for the ionized gas is defined as the ve- 
locity associated with the nucleus (region A), and it has 
a value of 3121±27 km s^^, consistent with Zhao et al. 
(1997) estimate based on the H92a radio recombination 
line (3120±50 km s~^). The ionized gas velocity field 
has a peak-to-peak variation of 290 km s~^. The neutral 
gas systemic velocity is 3057±34 km s^^, and the peak- 
to-peak velocity is about 216 km s~^. Although evident 
structural similarities are found between the ionized and 
neutral gas (i.e. same kinematics), in general the veloc- 
ity of the neutral gas presents a blueshift over the entire 
area with respect to the ionized gas (about -30±12 km 
s~^ over the studied regions). This systematic blueshift 
appears to be real taking into account the uncertainties. 
In that case, this is a marginal effect, with a lower value 
than the ones measured in ULIRGs using a slit (Mar- 
tin 2006; from about -60 km s^^ to -300 km s'^). In 
the present analysis the value obtained is derived from 
IFS data, i.e. has the advantage of the two dimensional 
distribution versus the use of the slit, where the measure- 
ments are limited to a restricted area, probably centered 
on the nuclear regions. 

Across the galaxy there are gradients of about 110 km 
s^^ kpc^^ in the ionized gas velocity (between regions lo- 
cated ^^2.3 kpc apart i.e. ~11"), higher than the typical 
values founded in ULIRGs (up to 60 km s~^ kpc~^ for an 
individual galaxy, Martin 2006) when using a slit of 20" 
long. In closer regions with large velocity variations the 
gradients can be 150 km s~^ kpc~^ (between regions I 
and northern of G, located about 1.4 kpc apart). In spite 
of these gradients, and as discussed below, there are no 
clear evidence for virialized simple rotating motions. 

The two-dimensional velocity dispersion map shows 
a different structure than the velocity field (Fig. 1181 
Right panel). For nucleus A, the a value is 108±11 km 
s~^, in well agreement with the derivation of Shier et 
al. (1996) based on the CO(2,0) near-infrared band at 
2.3 pLTn (135±21 km s~^). The maximum value for a 
(170±20 km s~^) is located in a region 0.28 kpc towards 
the north-east of A, where the [O i]A6300 emission peak 
is also located. Since [O i]A6300 trace shocks, this high a 
value is a hint of turbulence and shocks. The lowest ve- 
locity dispersion (lower limit 25±3 km s""'^) corresponds 
to the area dominated by H ii-like ionization, towards the 
east of the spiral arm (e.g. A7, see t l6.2|l . This indicates 
that the velocity dispersion is likely dominated by the in- 
ternal kinematics within the stellar clusters, rather than 
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by large regions projected along the line of sight. 

In this case it is not possible to state that there are 
clear indication of pure rotation, i.e. the photometric 
peak does not coincide with the velocity dispersion max- 
imum, and neither of them coincide with the possible 
kinematic center. The velocity maps show that even if 
pure rotation does not dominate the galaxy kinematics, 
rotational movements can be present. 

10.1.3. NGC 3690 

Like in the previous case, at a resolution of 0.18 kpc 
the structure of the ionized and neutral gas velocity fields 
obtained with the SB2 bundle are very similar (Fig. IT^ 
Left and Center panels). The largest receding velocities 
are detected towards the north-east of the field (C, C6), 
while the largest approaching velocities are located on B2 
and towards the south-west. The ionized gas systemic ve- 
locity is defined as the one of the nucleus Bl, with a value 
of 3040±27 km s^^. This value coincides with Zhao et 
al. (1997) measurements using the H92a radio recombi- 
nation line (3080±40 km s~^). For the ionized gas, the 
peak-to-peak velocity is 257 km s~^. The systemic ve- 
locity measured with the neutral gas is 2976±32km s~^, 
with a peak-to-peak variation of 217 km s~^. In both ion- 
ized and neutral gas, the velocity variations are similar 
to the ones measured in IC 694. As founded in ULIRGs 
(Martin 2006), there is a systemic blueshift of the neu- 
tral gas velocity with respect to the ionized gas (about 
-24±13km s^^ in the studied regions). Due to its value, 
lower than the results founded for ULIRGs (from about 
-60 km s"^ to -300 km s-\ Martin 2006), this blueshift 
is consider as a marginal effect. 

Across the galaxy there are gradients of about 70 km 
s^^ kpc^"'^ in the ionized gas (between C' and D3). This 
result is close to the gradients derived for ULIRGs using 
a slit (up to 60 km s~^ kpc~^ for an individual galaxy, 
Martin 2006). 

The two-dimensional structure of the velocity disper- 
sion ranges between 50±5 and 172±20 km s~^ (Fig. [T^ 
Right panel) ^ and present different characteristics than 
the two dimensional velocity fields. The nucleus Bl has 
the maximum velocity dispersion value of the entire sys- 
tem (172±20 km s'^), and regions with high velocity 
dispersion (~160 km s^^) are also detected about 1 kpc 
towards the north and north-west of Bl. Looking at the 
ionized gas velocity field, this secondary peaks of the ve- 
locity dispersion are located in regions where there are 
high velocity variations (from approaching to receding), 
and thus this a would represent a convolution of a large 
velocity gradient. Although the a-peak coincides spa- 
tially with the nucleus, and then it would be tracing the 
mass (large a implies that this is the most massive galaxy 
of the system), the presence of an AGN in Bl could pro- 
duce a broadening of the line width. 

Like in the case of IC 694 there are no clear evidences 
for ordered rotating motions; the photometric and ve- 
locity dispersion peaks do not coincide with the possible 
kinematic center. In conclusion, although the pure ro- 
tation does not dominate the kinematics of the galaxy, 
rotational process would be present. 

10.2. Where is the true nucleus of NGC 3690? 

The galaxy NGC 3690, shows variations on the surface 
brightness importance of the regions Bl and B2 with 



wavelength. While in the ultraviolet (Windhorst et al. 
2002) and in the optical (this work) B2 appears brightest 
than Bl, and therefore has been identified as the nucleus, 
the reverse is true in the near-IR (2.2 /im AAHOO) and 
mid-IR (Gallais et al. 2004). Besides, the optical kine- 
matical analysis shows that the a peak is located in Bl 
(see tilU.1.3|l . This fact along with the high extinction 
(Av=3.7), the hot dust component needed to explain 
the color excess, and the presence of an AGN, lead us 
to the conclusion that Bl is the true nucleus of the in- 
dividual galaxy NGC 3690, as previously discussed by 
Gallais et al. (2004) with mid-IR spectro-imaging, and 
concluded by Zezas et al. (2003) using Chandra X-ray 
observations. 

10.3. Dynamical mass derivation of the nuclei and 
brightest regions 

Until now the only mass derivation in this work has 
been done for the Ha-emitting young star forming re- 
gions. The central ionized gas velocity dispersion is a 
good tracer of the dynamical mass, thus it can be used 
as an alternative method for mass derivation. Assum- 
ing virialization, the dynamical mass of the nuclei and 
brightest regions can be derived as: 

Mdyn = 1.75 X lO^RhrnlH'lQ (3) 

where Rhm is the half mass radius in pc and a is the ve- 
locity dispersion in km s~^. Assuming a constant value 
of 1.75 (see Colina et al. 2005 for a detailed discussion), 
and that the effective radius (Re//) of the light distri- 
bution traces the half mass radius in each region, it is 
possible to infer the total dynamical mass. 

In the present treatment we assume that Ha is traced 
by Pack, that gives the real size of the emitting regions. 
In consequence the effective radius of the most impor- 
tant emitting sources is calculated over the HST Paa 
image; additional estimates using the FllOM image are 
done for comparison. There are some caveats related to 
the spatial resolution of the different set of data worth 
mentioning. For this mass estimate we are dealing with 
information at different spatial resolutions: on the one 
hand the velocity dispersion from the INTEGRAL SB2 
bundle data was measured at a 0'.'9 (~180 pc) resolu- 
tion. On the other hand the Re// over the most relevant 
Arp 299 regions (i.e. the ones resolved) was measured 
using the NICMOSl/FllOM image (O'.'043/pix, i.e. -9 
pc/pix) and NICM0S2/Paa image (O'.'076/pix, i.e. ~16 
pc/pix). The typical size of the Hli regions in LIRGs 
is 200-300 pc (Alonso-Herrero et al. 2002), i.e. several 
times larger than the effective radius of the stellar clus- 
ters and about the size of a SB2 fiber. Taking this into ac- 
count, the clusters will be point sources for INTEGRAL, 
and their a will be integrated on the fiber. 

For nucleus A we derive a dynamical mass of about 
1.2xlO^M0 (Table which is about 5 times smaller 
than the derivation done by Shier et al. (1996) using the 
velocity dispersion measured from the 2.3^m CO band; 
this discrepancy is due to a different estimate of the ra- 
dius (rco/rpoQ— 4). The calculated M^yn impHes that 
the young stellar mass content derived from L(IIa) in A 
is about 40%. The mean value for the dynamical mass of 
nucleus Bl is l.SxlO^M© (in good agreement with Shier 
et al. 1996), consequently the upper limit for the young 
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stellar mass content is about 15%. For B2 the mean 
value of the total mass is O.GxIO^Mq (the same as Shier 
et al. 1996). This would imply that about 100% of the 
mass in B2 is from young stellar population, indicating 
that the mass derivation done using L(Hq;) is probably 
overestimated. For the last two regions (CI and C), the 
Re// gives an upper limit for the dynamical mass (0.6 
and O.3xlO^M0 respectively, the latter according to the 
>0.1x10^Mq from AAHOO), that imphes young stellar 
mass contents of about 13% and 40% respectively. 

II. SUMMARY 

A detailed IFS study of the LIRG Arp 299 
(IC 694+NGC 3690) system covering a region of 7.0x6.2 
kpc has been presented. This work illustrates the flexi- 
bility and power of the IFS in studying in two dimensions 
and simultaneously the different physical properties, such 
as stellar populations, internal dust extinction, ioniza- 
tion structure and gas kinematics, of complex systems 
like Arp 299. In what follows only the main results of 
this study are summarized: 

1. The observed overall two dimensional stellar and 
ionized gas structure differs due to internal extinc- 
tion effects caused by the dust, and to the differ- 
ent spatial distribution of the diverse (AGN, stars, 
shocks) ionizing sources. The stellar continuum 
and emission line distribution peaks are not spa- 
tially coincident but separated up to 1.4 kpc. The 
peaks of the ionized gas traced by different emis- 
sion lines (e.g. hydrogen recombination, high ex- 
citation [OIII], and shock-tracer [01]) are also not 
coincident indicating that the ionizing mechanisms 
change on scales of about 200 pc or less. 

2. The region Bl is identified as the true nucleus 
of NGC 3690. The presence of an AGN in Bl 
is inferred from the detection of a high-excitation 
Seyfert-like conical structure of about 1.5 kpc in 
size, an opening angle of about 54 ° , and its apex 
being located in Bl. This region also presents the 
additional characteristics of a dust -enshrouded nu- 
cleus with high extinction (Av=3.7±0.6), red col- 
ors (F160W-F222M=2.06±0.1) indicating hot dust 
emission, and high velocity dispersion ((T=172±20 
km s~^). 

3. The observed near-IR colors reveal the presence of 
hot dust over the entire structure. According to he 
F160W-F222M color, the upper Hmit for the mean 
hot dust contribution is about 30% at 600 K for the 
majority of the star-forming regions. The nuclei A 
and Bl, and region CI are the most affected by the 



hot dust component, with an upper limit of about 
70%, 80% and 90% at 600 K respectively. 

4. The high surface brightness regions are mainly 
dominated by star formation, while the low surface 
brightness regions are dominated by LINER-like 
ionization process. The Ha-emitting star-forming 
regions present an extinction of between Ay = l— 6 
mag, with ages between 3.3 and 7.2 x 10^ years, and 
masses that range between about 6 and 680x10^ 
M0. The total contribution (extinction corrected) 
of these sources to the bolometric luminosity is 
about 45%. 

5. The structure of the ionized (Ha) and neutral 
(NaD) gas two dimensional velocity fields show 
similar features, indicating that both components 
are kinematically coupled on scales of hundreds of 
pc and larger. The velocity fields are very com- 
plex, and in general the galaxy does not appear 
to be dominated by ordered virialized simple (e.g. 
rotation) motions. The morphology of the veloc- 
ity field docs not correlate with the distribution of 
mass as traced by the continuum images, that is, 
no spatial coincidence between the kinematic and 
the photometric centers exists. 

On average, the interface region presents the lowest 
value of velocity dispersion (cr~45 km s~^), IC 694 
presents the intermediate values ((t~65 km s^^), 
and NGC 3690 the highest (<t~120 km s'^). This 
indicate that NGC 3690 is the warmest body of 
the entire system, and probably the most massive. 
On the contrary, the low velocity amplitude and 
velocity dispersion of the interface suggests that the 
ionized gas is slowly rotating or in a close quiescent 
phase. 

The velocity dispersion measurements over the 
nuclei and the brightest emitting regions, imply 
a mean value for the dynamical mass of about 
1.2x10^ Mq and 1.5 xlO^ for nuclei A and Bl re- 
spectively, within a mean radius of 58 pc and 28 
pc. For regions B2, CI and C the mass ranges be- 
tween -0.2x10^ M© and ~0.6xl0^ M©, within a 
radius of about 30 pc. 
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TABLE 1 

Properties of the regions of the IC 694+NGC 3690 interface 



Reg 


Dist'' 
(kpc) 


(circsec) 


(circsec) 




Fobs (Ha)'' 


(mag) 




log([OIII]/H/3)d 


log([OI]/Ha)'= 


log([NII]/Ha)' 


log([SII]/Ha)f5 


Kl 


2.97 


-17.6 


-12.7 


1.7 


2.2 


0.27 


-1.07 


-0.34 


-0.37 


K2 


3.81 


-21.6 


-13.7 


1.1 


0.9 


0.07 


-1.04 


-0.36 


-0.38 


K3 


1.28 


-9.6 


-10.7 


0.6 


2.0 


0.22 


-0.59 


-0.21 


-0.25 


K4 


2.49 


-13.0 


-1.8 


2.7 


1.6 


0.04 


-0.83 


-0.41 


-0.31 


K5 


3.53 


-14.5 


3.9 


10.7 


2.9 


0.14 


-1.44 


-0.54 


-0.64 


K6 


3.71 


-13.1 


6.0 


5.7 


2.3 


-0.05 


-1.27 


-0.48 


-0.54 


K7 


2.72 


-10.6 


1.9 


4.0 


2.0 


0.01 


-1.10 


-0.47 


-0.45 


K8 


4.07 


-9.0 


9.8 


1.4 


2.7 


0.14 


-0.71 


-0.32 


-0.23 


K9 


3.52 


-2.6 


7.8 


3.3 


1.8 


-0.01 


-0.92 


-0.37 


-0.36 


KIO 


5.26 


-21.5 


-9.0 


4.2 


1.0 


0.17 


-1.36 


-0.53 


-0.56 


Kll 


1.44 


-9.1 


-5.1 


1.6 


1.5 


0.02 


-0.82 


-0.39 


-0.31 


K12 


3.18 


-16.8 


-17.4 


7.6 


2.0 


0.13 


-0.98 


-0.35 


-0.37 



Note. — (a): Distance relative to the nucleus A. (b): Flux measured (xlO^^^crg s^'cm^^) using a circular aperture of 3" Uncertainties 
10%-15%. (c): Internal extinction derived from the Ha/H/3 ratio. Average error in Ay is 0.7 magnitudes, (d): Average error ±0.11 (e): 
Average error ±0.06 (f): Average error ±0.05 (g); Average error ±0.06. In all cases, emission line ratios are derived using extinction 
corrected flux values. 



14 



M. Garcia-Marin et al. 



TABLE 2 

Properties of the regions of IC 694 



Reg 


Dist^' 
(kpc) 


(circsec) 


(circsec) 


Foi„(Ha)b 


(mag) 


log([OIII]/H/3)d 


log([OI]/Ha)<= 


log([NII]/Ha)f 


log([SII]/HQ)f5 


A 


0.00 


0.0 


0.0 


2.9 


2.9 


-0.14 


-1.03 


-0.31 


-0.41 


Al 


0.46 


1.5 


-1.6 


2.6 


2.3 


-0.22 


-1.29 


-0.40 


-0.50 


A2 


0.54 


-1.0 


-2.4 


3.1 


2.0 


-0.14 


-1.51 


-0.52 


-0.60 


A3 


1.37 


3.3 


-5.7 


2.3 


1.9 


-0.03 


-1.61 


-0.53 


-0.59 


A4 


1.51 


3.8 


-6.3 


3.0 


2.1 


-0.02 


-1.73 


-0.55 


-0.66 


A5 


1.06 


-2.1 


-4.7 


6.0 


1.9 


0.10 


-1.56 


-0.54 


-0.62 


A6 


1.00 


-3.6 


-3.9 


4.0 


2.1 


0.04 


-1.46 


-0.52 


-0.56 


A7 


0.99 


-2.8 


-3.9 


6.2 


2.3 


0.09 


-1.61 


-0.56 


-0.66 


F 


1.58 


6.5 


-4.1 


6.7 


3.2 


0.10 


-2.06 


-0.52 


-0.92 


G 


1.10 


1.1 


-5.2 


6.2 


1.9 


0.05 


-1.84 


-0.54 


-0.74 


HI 


0.48 


2.3 


-0.3 


6.4 


2.8 


-0.15 


-1.41 


-0.47 


-0.63 


H2 


0.34 


1.6 


0.5 


4.7 


3.4 


-0.06 


-0.99 


-0.32 


-0.46 


I 


0.42 


-1.8 


-0.9 


6.0 


1.6 


-0.09 


-1.60 


-0.52 


-0.65 


T1 


1 45 


-70 


8 


0.8 


9 4 


-mn 


-^ 9.H 


-0.43 


-n /I?! 



Note. — (a): Distance relative to the nucleus A. (b): Flux measured (xlO^^'^erg s~^cm~^) using a circular aperture of l" Uncertainties 
10% - 15%. (c); Internal extinction derived from the Ha/H/3 ratio. Error in Ay is 0.5 magnitudes, (d): Average error ±0.08 (e): Average 
error ±0.04 (f): Average error ±0.05 (g): Average error ±0.04. In all cases, emission line ratios are derived using extinction corrected flux 
values. 
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TABLE 3 

Properties of the regions of NGC 3690 



Reg 


Dist^ 
(kpc) 


(circsec) 


A(5^ 
(circsec) 




(mag) 


log([OIII]/H/3)'' 


log([OI]/Ha)<= 


1 /riVTXXl /XT \"f 

log([NII]/Ho:)' 


log([SII]/Ha)^ 


Bl 


0.00 


0.0 


0.0 


9.7 


3.7 


0.55 


-1.07 


-0.29 


-0.54 


B2 


0.49 


-1.5 


1.9 


4.0 


3.9 


0.43 


-1.37 


-0.35 


-0.74 


Bll 


0.22 


0.27 


-1.0 


5.7 


2.0 


0.10 


-1.37 


-0.41 


-0.52 


B16 


0.26 


-0.5 


-1.1 


10.6 


1.1 


-0.02 


-1.68 


-0.47 


-0.71 


B21 


0.82 


-1.0 


3.9 


3.1 


2.6 


0.02 


-1.54 


-0.48 


-0.77 


B22 


0.63 


-2.9 


1.0 


3.2 


2.1 


0.02 


-1.51 


-0.41 


-0.68 


B25 


0.85 


-3.9 


1.3 


7.2 


1.7 


-0.02 


-1.88 


-0.47 


-0.80 


CI 


1.83 


-2.4 


8.6 


28.1 


3.0 


0.13 


-1.87 


-0.49 


-0.96 


C2 


1.80 


-3.0 


8.2 


33.4 


2.1 


0.11 


-1.92 


-0.46 


-0.94 


C3 


1.75 


-3.9 


7.6 


16.7 


1.0 


0.01 


-1.80 


-0.44 


-0.78 


C4 


1.39 


-4.2 


5.3 


5.7 


0.5 


-0.03 


-1.80 


-0.47 


-0.67 


06 


1.55 


1.0 


7.5 


5.2 


1.9 


0.02 


-1.74 


-0.53 


-0.70 


c' 


1.96 


2.7 


9.5 


0.9 


3.8 


0.06 


-1.40 


-0.45 


-0.57 


D3 


1.49 


-7.2 


-0.9 


1.3 


1.3 


-0.03 


-1.60 


-0.44 


-0.55 


D4 


1.13 


-5.5 


0.1 


2.6 


2.0 


-0.05 


-1.44 


-0.41 


-0.57 



Note. — (a): Distance relative to the nucleus Bl. (b): Flux measured (xlO^^'^erg s~^cm~'^) using a circular aperture of 0''9. Uncer- 
tainties 10% - 15% (~25% for B2).(c): Internal extinction derived from the Ha/H/3 ratio. Error in Ay is 0.5 magnitudes (0.6 for Bl and 
C , and ~ 0.7 for B2). (d): Average error ±0.13 (e): Average error ±0.13 (f): Average error ±0.05 (g): Average error ±0.05. In all cases, 
emission line ratios are derived using extinction corrected flux values. 
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TABLE 4 

PAa/Ha-BASED INTERNAL EXTINCTION ESTIMATES 



Region 


f(Paa) 
(xlO~l^erg cm~2s 


f(Paa)/f(HQ) 


(mag) 


A 


167 


5.8 


6.1 


Al 


2.74 


0.1 




A2 


5.68 


0.2 


0.9 


A3 


4.37 


0.2 


0.9 


A4 


9.62 


0.3 


1.5 


A5 


8.41 


0.1 




A6 


2.73 


0.1 




A7 


10.5 


0.2 


0.9 


HI 


12.6 


0.2 


0.9 


H2 


17.3 


0.4 


2.0 


I 


5.11 


0.1 




11 








Bl 


66.7 


0.7 


2.8 


B2 


6.10 


0.2 


0.9 


Bll 


9.77 


0.2 


0.9 




±0.0 




n Q 


B21 


4.14 


0.1 




B22 


5.15 


0.2 


0.9 


B25 


11.0 


0.2 


0.9 


CI 


85.3 


0.3 


1.5 


C2 


76.3 


0.2 


0.9 


C3 


10.0 


0.1 




C4 


4.70 


0.1 




C6 


9.16 


0.2 


0.9 


C 


38.1 


4.1 


5.6 


D3 


0.15 


0.01 




D4 


4.81 


0.2 


0.9 


Note. — Fluxes of Pao measured over the selected regions using the HST (F190N-F187N) image, 
magnitudes. Some values axe discarded due to the noise alTecting Paa and Ha low-surface brightness. 


Average uncertainty in Ky is 0.7 
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TABLE 5 
Observed magnitudes and colors 



Region 


niF814W 




F110M-F160W 


F160W-F222M 


F814W-F110M 




16.85 


13.72 


1.33 


1.36 


1.80 


Al 


18.53 


17.59 


1.07 


1.04 


-0.13 


A2 


-LO.O i 




n Qi 


n f^7 


u.ou 


A3 


18.40 


17.75 


0.81 


0.65 


-0.15 


A4 




1 fi no 


1 




-U.O J. 


A5 


18.19 


17.61 


0.72 


0.63 


-0.14 


A6 


17.87 


17.85 


0.70 


0.54 


-0.68 


Bl 


15.67 


13.93 


1.23 


2.06 


0.51 


Bll 


17.30 


16.31 


0.88 


0.64 


0.11 


B16 


16.90 


18.12 


0.19 




-1.42 


B2 


14.09 


13.19 


0.70 


0.69 


0.20 


B21 


16.16 


15.35 


0.79 


0.55 


0.03 


B22 


17.78 


16.49 


0.93 


0.75 


0.36 


B25 


18.01 


17.67 


0.94 


0.62 


-0.56 


CI 


17.04 


14.68 


1.71 


2.67 


0.66 


C2 


16.63 


16.57 


0.67 




-0.61 


C3 


17.04 


16.37 


0.70 


0.52 


-0.02 


C4 


17.24 


16.71 


0.82 


0.44 


-0.29 


D3 


18.74 


18.37 


0.70 


0.46 


-0.33 



Note. — 



Errors in the colors are ±0.10 and ±0.14 for the bright and faint sources respectively. 
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TABLE 6 

Properties of the stellar populations in regions of IC 694 and NGC 3690 



Region 


L(Ha) 


TJITI 7/TT ^ \ 


Age 




Mass(Mpieow) 




(xlO*"ergs ^) 


(A) 


(xlO® years) 




(xIO^Mq) 






201 


5.3 






Al 


4.8 


217 


5.2 


17.6 


6.0 


A2 


4.6 


292 


5.0 


13.0 


6.4 


A3 


3.2 


297 


5.0 


9.0 


4.8 


A4 


4.7 


382 


4.8 


10.0 


4.0 


A5 


7.3 


307 


5.0 


21.0 


5.3 


A6 


6.1 


197 


5.3 


26.8 


5.0 


A7 


11.3 


277 


5.1 


36.5 




p 


28.6 


760 


3.3 


21.8 




G 


7.9 


654 


3.4 


6.2 




HI 


18.6 


416 


4.7 


34.0 




H2 


22.5 


314 


4.9 


57.0 




J 


5.8 


220 


5.2 


21.5 




IX 


1.6 


86 


5.9 


12.9 




Bl 


62.3* 


240* 


5.2* 


231.0* 


200.0* 


B2 


29.8 


22 


7.2 


678.0 


300.0 


Bll 


7.9 


200 


5.3 


34.6 


14.0 


B16 


7.1 


311 


5.0 


20.0 


3.1 


B21 


7.6 


62 


6.1 


74.0 


58.0 


B22 


4.8 


138 


5.6 


27.3 


10.0 


B25 


7.7 


449 


4.6 


13.7 


5.6 


CI 


93.8 


680 


3.4 


74.0 


81.0 


C2 


53.4 


645 


3.5 


43.8 


9.8 


C3 


9.7 


440 


4.7 


17.7 


9.8 


C4 


2.1 


156 


5.4 


10.3 


8.4 


C6 


6.8 


609 


3.6 


5.7 




c' 


5.2 (31.0) 


245 


5.2 


19.0 (115.0) 




D3 


1.0 


116 


5.7 


6.4 


3.2 


D4 


3.6 


360 


4.8 


7.8 





Note. — (a) Luminosity corrected for extinction using Ha /Hp ratio; the extinction using Paa/Ha ratio is given for regions A and C 
in brackets, (b) Mass derived from the Ha: luminosity corrected by extinction using Ha/Hj8 ratio, and assuming the age listed in column 
4. In brackets estimate using the Paxx/Ha extinction-corrected data, (c) Mass derived from the absolute magnitude yipieow corrected by 
extinction using Ha/H^S ratio, and assuming age listed in column 4. In bra<;kets estimate using the Paa/Ha extinction-corrected data. * 
means AGN contribution. 
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TABLE 7 
Kinematic of the interface region 



Region 


V(Ha)^ 


V(Nai)'' 


<7(Ha)<= 




(km s-l) 


(km s-l) 


(km s-l) 


Kl 


-56 


-29: 


61 


K2 


-32 




58 


K3 


-65 


30 


80 


K4 


13 


31 


68 


K5 


54 


101 


35 


K6 


54 


115 


20 


K7 


36 


78 


30 


K8 


84 


42 


70 


K9 


20 


-6: 


65 


KIO 


54 


86: 


100 


Kll 


-2 


26 


19 


K12 


-78 




36 



Note. — The mean errors for the velocity measurements in Ha and Nai for this bundle are 27 km s ^ and 60 km s ^ respec- 
tively.":" means large uncertainties. Velocities referred to nucleus (region A) of galaxy IC 694. (a) VA(Ha)=3121±27kms~^.(b) 
VA(Nai)=3057±34kms-i. (c) Relative error in a is ~13% 
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TABLE 8 
Kinematic of IC 694 regions 



Region 


V(Ha)^ 


V(Nai)'' 


<7(Ha)<= 




(km s-l) 


(km s-l) 


(km s-i) 


A 








108 


Al 


-86 


-47 


73 


A2 


-7 


22 


64 


A3 


-134 


-86 


36 


A4 


-151 


-96 


50 


A5 


-67 


-25 


45 


A6 


-44 


-37 


40 


A7 


-57 


-26 


41 


F 


-163 




72 


G 


-121 


-77 


43 


HI 


-130 


-89 


61 


H2 


-117 


-93 


109 


T 


fin 


Q8 


fil 


11 Note. — 


Velocities referred to nucleus (re|;Ton A), (a) V>t(Ha)=3121±27kms 


;- 1 . (b) (Na i) =3057±34 km s" ^ . 


(c) Relative ISror in 


u is ~10% 
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TABLE 9 
Kinematic of NGC 3690 regions 



Region 


V(Ha)^ 
(km s~l) 


V(Nal) 
(km s ^) 


u(tia) 
(km s~i) 


Bl 








172 


B2 


-116 


-55 


110 


Bll 


38 


50 


143 


B16 


-46 


-15 


117 


B21 


-18 


9 


128 


B22 


-121 


-48 


78 


B25 


-112 


-64 


54 


CI 


32 


65 


93 


C2 


16 


51 


99 


C3 


13 


36 


90 


C4 


-24 


8 


94 


C6 


115 




82 


c' 


101 




69 


D3 


-101 


-40 


67 


D4 


-112 


-52 


54 


Note. — Velocities referred to nucleus (regions Bl). (a) Vsi(Ha)=3040±27kms 
in a is ~10% 


i-^(b) Vsi(Nal)=2976±32kms 


. (c) Relative error 
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TABLE 10 

Dynamical mass estimates 



Region 


Re//" 






(kpc) 


(xIOSMq) 


A 


0.054-0.062 


1.10-1.3 


Bl 


0.035-0.022 


1.81-1.1 


B2 


0.030-0.025 


0.63-0.5 


CI 


0.022-0.042 


0.33-0.6 


c' 


0.023-0.043 


0.19-0.3 



Note. — Dynamical mass estimates for the nuclei of the system (A and Bl) and main Hcj-selected emitting regions {B2, CI and C'). 
The stellar and gas effective radius have been derived for the HST FllOM (first column) and for F190N (Paa) (second column). In both 
cases the dynamical mass has been derived using the Ho velocity dispersion values. See H10.3l for details. 
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Fig. 1. — Archival HST image of the interacting system Arp 299 obtained with the WFPC2 F814W filter. The rectangular overlays 
indicate the effective area covered by the three pointings with the two INTEGRAL configurations used (i.e SB2 small rectangles, SB3 
large rectangle). We have marked the positions of the regions under study (see § 3 and Tables 1, 2, and 3). The location of the recently 
discovered near-IR SN is shown too (Mattila et al. 2005). The image is shown on a logarithmic scale. [See the electronic edition of the 
Journal for a color version of this figure.] 
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Fig. 2.— Average spectra of the Interface Region (Kl to K12), IC 694 (A, HI and H2) and NGC 3690 (Bl and Bll). The spectral range 
shown is 4500-6900A. See Fig. [Tland Tables [Tl 1^ and O for location of the individual emitting sources. 
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Fig. 3. — Emission-line free stellar continuum and emission line images for Arp 299 obtained with INTEGRAL SB3 bundle. NGC 3690 
(Bl and B2) is located to the west and IC 694 (A) to the east. The contours represent the red continuum at A6460. The HST /WFPC2 
F814W image is shown for comparison, where several regions of interest discussed in the text arc marked. All the images are shown on a 
logarithmic scale, and the color code represents different intensity levels in each of the individual maps. Orientation is north up, east to 
the left. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 4. — Emission-line free stellar continuum and emission line images for IC 694 obtained with INTEGRAL SB2 bundle. The contours 
represent the red continuum at A6460. The HST/WFPC2 F814W image is also shown for comparison, with several regions of interest 
marked on it. All the images are shown on a logarithmic scale, and the color code represents different intensity levels in each of the 
individual maps. Orientation is north up, east to the left. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 5. — Emission-line free stellar continuum and emission line images for NGC 3690 obtained with INTEGRAL SB2 bundle. The 
contours represent the red continuum at A6460. The HST/WFPC2 F814W image is also shown for comparison, with several regions of 
interest marked on it. All the images are shown on a logarithmic scale, and the color code represents different intensity levels in each of 
the individual maps. Orientation is north up, east to the left. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 6. — Left panel: Extinction corrected H« map of the Arp 299 system (including the interface region), shown on a logarithmic scale. 
Right panel: Extinction map of the interface region obtained from the line ratio Ha/H/3, shown on a linear scale. Contours are as in Fig. 
131 [See the electronic edition of the Journal for a color version of this figure.] 




Fig. 7.— IC 694. Top left panel: HST/NICMOS continuum subtracted Pa« emission (F190N-F187N). Top right panel: HST/NICMOS 
H-band continuum. Bottom left panel: INTEGRAL Ha extinction corrected map. Bottom right panel: Extinction map obtained with 
Ha/H/3 line ratio. As a reference several regions in the two fields are marked, and the red continuum contours are superimposed on the 
interpolated maps. The NICMOS images are shown covering the same FoV as the INTEGRAL data. All the images are shown on a 
logarithmic scale except for the extinction map, where a linear scale has been applied. [See the electronic edition of the Journal for a color 
version of this figure.] 
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Fig. 8.— NGC3690. Top left panel: HST/NICMOS continuum subtracted Paa emission (F190N-F187N). Top right panel: HST/NICMOS 
H-band continuum. Bottom, left panel: INTEGRAL Ha extinction corrected map. Bottom right panel: Extinction map obtained with 
Ha/H/3 line ratio. As a reference several regions are marked, and the red continuum contours are superimposed on the interpolated maps. 
The NICMOS images are shown covering the same FoV as the INTEGRAL data. All the images are shown on a logarithmic scale except 
for the extinction map, where a linear scale has been applied. [See the electronic edition of the Journal for a color version of this figure.] 




Fig. 9. — Optical emission line diagnostic diagrams for the studied regions of Arp 299. Filled symbols represent nuclei A and Bl, and 
region B2. Note that Bl, as observed with the SB2 bundle, is classified as Seyfert according to these diagrams. The open symbols are the 
regions given in Tables 1, 2, and 3. Typical uncertainties in the line ratios of each individual galaxy are shown. [See the electronic edition 
of the Journal for a color version of this figure.] 
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Fig. 10. — Diagnostic diagram and line ratio maps derived from INTEGRAL/SB3 data. Upper panel: Diagnostic diagram with all the 
data from the SB3 bundle displayed. The symbol code is as follows: triangles represent values within the ionization cone, and at a distance 
of less than 2 kpc from Bl. Stars represent the zone within the cone at distances of 2 to 4 kpc from Bl. Crosses mark the rest of the 
regions in the interface. Filled circles represent regions in IC 694 and NGC 3690 (defined as the regions within the iso-contours). Bl has 
been included in the ionization cone, and labeled as a triangle. Center panel: [O in] A5007/H/3 map. Nuclei A and Bl are marked. The 
white lines indicate the ionization cone. The transversal line to the axis of the cone separates the inner (with Seyfert-like regions) and outer 
(LINER-like) sections of the cone. Lower panel: [Ol]A6300/Ha map. The labels and overlays as in the center panel. [See the electronic 
edition of the Journal for a color version of this figure.] 
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Fig. 11. — Emission line diagnostic diagram and excitation maps based on line ratios for IC 694. Upper panel: Diagnostic diagram 
obtained for the individual SB2 spectra. Center and Lower panels: Diagnostic maps derived from reconstructed emission line maps. As a 
reference we have marked several regions and superimposed red continuum contours. [See the electronic edition of the Journal for a color 
version of this figure.] 
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Fig. 12. — Emission line diagnostic diagram and excitation maps based on line ratios for NGC 3690. Upper panel: Diagnostic diagram 
obtained from the individual SB2 spectra. Center and Lower panels: Diagnostic maps derived from reconstructed emission line maps. As 
a reference we have marked several regions and superimposed red continuum contours. [See the electronic edition of the Journal for a color 
version of this figure.] 
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Fig. 13. — Relation between velocity dispersion and [Ol]A6300/Ha ratio for the three INTEGRAL pointings, i.e. interface region (Left), 
IC 694 {Center panel) and NGC 3690 (Right panel). In all cases the relation derived by Monreal-Ibero et al. (2006) for a set of ULIRGs 
is shown (solid line). The shaded area indicate the region covered by extra- nuclear regions of ULIRGs (Monreal-Ibero et al. 2006). The 
horizontal dashed line marks the average values separating Hll and LINER-like ionization. Left panel: The symbol code as in Fig. 10. 
Center and Right panels: Filled symbols indicate the fibers of the studied individual regions listed in Tables 121 and l3l {See the electronic 
edition of the Journal for a color version of this figure] 
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Fl10M-FieOW FieOW-F222M 

Fig. 15. — Plots representing the results of the stellar evolutionary models (SB99, Geneva High Isochrones, Salpeter IMF, using upper and 
lower mass boundaries of 0.1 and 100 M0) compared with the measured values. The long tracks represent the evolution of instantaneous 
bursts as a function of age (from top to bottom, black crosses mark 4, 6, 8 and 10 Myr). Short tracks correspond to continuous bursts using 
a SFR of 1 Mq year~^. In both cases from left to right the extinction ranges from Ay=0 to Av=6, in steps of 2 magnitudes. The filled 
circles are the measured values, and the stars represent the extinction corrected values from the Ha/H/3 ratio. Left panel: F110M-F160W 
color vs Ha equivalent width. Right panel: F160W-F222M color vs Ha equivalent width. The long horizontal lines represent the colors 
expected from a extinction-free stellar population of 6x10^ years with an increasing K-band flux contribution of hot dust at different 
temperatures (from 0% to 100%, increasing 10% on each step). In both cases nucleus Bl is included for completness. [See the electronic 
edition of the Journal for a color version of this figure.] 




Fig. 16. — Results of the stellar evolutionary models (SB99, Geneva High Isochrones, Salpeter IMF, using upper and lower mass 
boundaries of 0.1 and 100 Mq) compared with measured values. Long tracks are the output for instantaneous bursts as a function of age 
(from top to bottom, black crosses mark 4, 6, 8 and 10 Myr), and different cluster mass. The filled circles represents measured absolute 
F160W magnitude, and the stars represent the extinction corrected values from the Ha/H/3 ratio. [See the electronic edition of the Journal 
for a color version of this figure.] 
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Fig. 17. — Arp 299 interface region kinematics. Velocity field of the ionized gas as traced by the Ha line {Left panel), neutral gas as 
traced by the Na I line {Center panel) and velocity dispersion of the ionized gas {Right panel). The velocities are referred to the nucleus A 
as observed with the higher spatial resoluti on S B2 bundle, i.e. 3121±27 km s~^ for the ionized and 3057it34 km s~^ for the neutral gas. 
Lines depicted the ionization cone as in Fig llUI [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 18. — IC 694 kinematics. Velocity field of the ionized gas as traced by the Ha line {Left panel), neutral gas as traced by the Na 
I line {Center panel) and velocity dispersion of the ionized gas {Right panel). The velocities are referred to the nucleus A, for which the 
ionized gas velocity is 3121±27 km s~^ and the neutral gas velocity is 3057±34 km s~^. [See the electronic edition of the Journal for a 
color version of this figure.] 
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Fig. 19. — NGC 3690 kinematics. Velocity field of the ionized gas as traced by the Ha line {Left panel), neutral gas as traced by the Na 
I line {Center panel) and velocity dispersion of the ionized gas {Right panel). The velocities are referred to the nucleus Bl, for which the 
ionized gas velocity is 3040it27 km s~^ and the neutral gas velocity is 2976±32 km s~^. Note that the nucleus Bl presents the the highest 
cr value. [See the electronic edition of the Journal for a color version of this figure.] 



